Evidence of a new narrow resonance decaying to Xci7 in B ^ XcilK 
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We report measurements of B — >■ xa'yK and Xc2'yK decays using 772 x 10*' BE events collected 
at the T(4S') resonance with the Belle detector at the KEKB asymmetric-energy e^e^ collider. 
Evidence of a new resonance in the Xci7 final state is found with a statistical significance of 3.8 a. 
This state has a mass of 3823.1 ± 1.8(stat) ± 0.7(syst) MeV/c^, a value that is consistent with 
theoretical expectations for the previously unseen l'^D2 cc meson. We find no other narrow resonance 
and set upper limits on the branching fractions of the X(3872) — >■ Xal s.nd Xc27 decays. 

PACS numbers: 13.25.Hw, 13.20.Gd, 14.40. Pq 



During the last decade, a number of new charmonium 
(cc)-like states were observed, many of which are candi- 
dates for exotic states [H. The first of these, the X(3872), 
has been observed by six different experiments in the 
same final state 0-0| • A recent update from Belle Q and 
LHCb results in a world average mass at 3871.68±0.17 
MeV/c2 d and a stringent upper bound on its width 
(r <_1.2 MeV) The proximity of its mass to the 
£)*0£)0 threshold makes it a good candidate for a DD* 
molecule [13] ■ Other, but alternative models have been 
proposed, such as a tetraquark or a ccg hybrid me- 
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Radiative decays can illuminate clearly the nature of 
hadrons. For example, the observation of X(3872) — > 
J/ipj confirmed the C-even parity assignment for the 



X(3872) [I3j, OJ- The similar decay mode X(3872) 



ip'^ could help to distinguish whether the X(3872) is 
charmoniuin, a molecular state or a mixture of the 
two Belle reported no significant signal in 

X(3872) — >■ ' 0^7 , in disagreement with positive evidence 
from BaBar The X(3872) Xci7 and Xc27 decays 
arc forbidden by C-parity conservation in electromag- 
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netic processes. However, if the X(3872) is a tetraquark 
or a molecular state, it may have a C-odd partner, which 
could decay into Xci7 and Xc27 final states |17l. Il8|. 

In the charmonium family, the observation of a U-wave 
cc meson and its decay modes would test phenomenolog- 
ical models [l^ E^l- The as- yet undiscovered 1^D2 cc 
(■02) and l^Ds cc (tps) states are expected to have sig- 
nificant branching fractions to XciJ and Xc27, respec- 



tively [21|, 1221. D-wave cc states and their properties 



were predicted long ago but remain unconfirmed [2l|, [2^ . 
The E705 experiment reported an indication of a 1^D2 
state in tt^N J/?^+7r^+ anything, with a mass of 
3836 ± 13 MeV/c^ [23|; however, the statistical signifi- 
cance of this result was below the threshold for evidence. 

In this letter, we report measurements of _B — Xci7^ 
and B -^_Xc2lK decays, where the Xci and Xc2 decay 
to J/V'7 [2J1. These results are obtained from a data 
sample of 772 x 10^ BB events collected with the Belle 
detector [2^ at the KEKB asymmetric-energy e+e~ col- 
lider operating at the T(4S') resonance [2^. The Belle 
detector is a large-solid-angle spectrometer, which in- 
cludes a silicon vertex detector (SVD), a 50-layer cen- 
tral drift chamber (CDC), an array of aerogel thresh- 
old Cherenkov counters (ACC), time-of- flight scintilla- 
tion counters (TOF), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux return located outside the coil 
is instrumented to detect K'j^ mesons and identify muons 
(KLM). 

The J /ij] meson is reconstructed via its decays to 
(I = e or /z). To reduce the radiative tail in the e'^e~ 
mode, the four-momenta of all photons within 50 mrad 
with respect to the original direction of the e"*" or e~ 
tracks are included in the invariant mass calculation, 
hereinafter denoted as Me+e-('y). The reconstructed in- 
variant mass of the J/'0 candidates is required to sat- 



Charged kaons arc identified by combining informa- 
tion from the CDC, TOF, and ACC systems. The kaon 
identification efficiency is 89% while the probability of 
misidentifying a pion as a kaon is 10%. mesons are 
reconstructed by combining two oppositely charged pi- 
ous with the invariant mass Af^r+ir- lying between 482 
MeV/c^ and 514 MeV/c^. The selected candidates are re- 
quired to satisfy the quality criteria described in Ref. . 

To reconstruct B candidates, each Xcj |2§| is com- 
bined with a kaon candidate and a photon having > 
100 MeV (and not used in the reconstruction of Xcj)- 
If the invariant mass of any photon pair that includes 
this photon is found to be consistent with a 7r° (i.e., 117 
MeV/c^ < M^^ < 153 MeV/c^), this photon is rejected. 
Among the events containing at least one Xcj candidate, 
9.0% have multiple Xcj candidates. In such cases, the 7 
forming the YcJ candidate with mass closest to the Xci or 
Xc2 masses [9| is not used as the additional photon. This 
treatment suppresses reflections from the Xcj daughter 
photons. 

The B candidate is identified by two kinematic 
variables: the beam-constrained mass (A/bc = 
and the energy difference {AE 



beam 



El 



-^boam)- Here, i?bcam i^ ^hc run-depcndeut beam energy, 
and E*Q and p*^ are the reconstructed energy and mo- 
mentum, respectively, of the B meson candidates in the 
center-of-mass (CM) frame. Candidates within a IS.E 
window of [-28, 30] MeV and with Mbc > 5.23 GeV/c^ 
are selected. Of these, 9.8% (6.4%) have multiple candi- 
dates in the B^ — s- XcilK^ [B^ Xc2jK^) mode; we 
select the B candidate with AE closest to zero. In order 
to improve the resolution in M^^^^, we scale the energy 
of the 7 so that AE is equal to zero. This corrects for 
incomplete energy measurement in the ECL. To suppress 
continuum background, events having a ratio R2 of the 
second to zeroth Fox- Wolfram moments [s^l above 0.5 
arc rejected. 



isfy 2.95 GeV/c^ < M^. 



-(7) 



< 3.13 GeV/c^ or 3.03 The Af^^,^ and Afbc projections for the B± ^ XalK^ 



GeV/c^ < A/^+^- < 3.13 GeV/c^. For the selected J/i^ 
candidates, a vertex-constrained fit is applied and then 
a mass-constrained fit is performed in order to improve 
the momentum resolution. The Xci and Xc2 candidates 
are reconstructed by combining J /if) candidates with a 
photon having energy {E^) larger than 200 MeV in the 
laboratory frame. Photons are reconstructed from the 
energy deposition in the ECL by requiring no matching 
with any extrapolated charged track. To reduce the back- 
ground from 7r° — 7^ 77, we use a likelihood function that 
distinguishes an isolated photon from 7r° decays using 
the photon pair invariant mass, photon laboratory en- 
ergy and polar angle |27[. We reject both 7's in the pair 



signal candidates are shown in Fig.[Tl where a 7/;' — Xci7 
signal is evident. In addition, there is a significant narrow 
peak at 3823 MeV/c^, denoted hereinafter as X(3823). 
No signal of ^(3872) — )• Xci7 is seen. We extract the 
signal yield from a two-dimensional unbinned extended 
maximum likelihood (2D UML) fit to the variables M^^^^ 
and Afbc- The likelihood function is 



N 



4=1 k 



if the 7r° likelihood probability is larger than 0.7. The 
reconstructed invariant mass of the Xci (Xc2) is required 
to satisfy 3.467 GeV/c^ < Mj/^^ < 3.535 GeV/c^ (3.535 
GeV/c^ < Mj/^^ < 3.611 GeV/c^). A mass-constrained 
fit is applied to the selected Xci and Xc2 candidates. 



where k is the index for the signal and background con- 
tributions, Uk and Pk are the yield and probability den- 
sity function (PDF), respectively, for the fc-th compo- 
nent, mxi and A^bi are the Af^^,^ and Afbc values in the 
i-th event, respectively, and N is the total number of 
candidate events in the data. 

The resolution in Af^^j^ (Afbc) is parameterized by a 
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sum of two Gaussians (Gaussian and logarithmic Gaus- 
sian |3l|). MC studies show that the resolutions in both 



-^x<!j7 ^^"i -^bc for a narrow resonance in the mass range 
3.8 GeV/c^ < M^^j^ < 4.0 GeV/c^ are in good agree- 
ment with those for The parameters of the resolution 
functions are determined from the MC simulation that is 
calibrated using the B"^ — >• tp'{^ Xcil)K'^ signal. We 
take into account the ip' natural width |9j by convolving 
the Breit-Wigner function and the resolution function; 
for the X{3823) and X(3872), zero natural widths are 
assumed. The two-dimensional PDF is a product of the 
one-dimensional distributions. 

For — > ip'{^ Xcil)K^ decay, the mean and width 
of the core Gaussian are floated and the remaining pa- 
rameters are fixed according to MC. To fit the B^ 
X(3823)(— ^ Xcij)K^ signal, we float the mean of the 
core Gaussian but constrain the detector resolution by 
using the ip' signal results after taking into account the 
difference estimated from the signal MC study. For Mbc 
the parameters are fixed to those found for the ip' , in ac- 
cordance with expectations based on the MC simulation. 
We constrain the mass difference (from Ref. [§]) and the 
sigma difference (from MC simulations) between ^' and 
X(3872). To fit B° ip'i^ Xcil)K°, B° X(3823)(^ 
Xcil)K°, B ^ ^'(^ Xc2j)K, B ^ X(3823)(-> Xc2l)K 
and B — >• X(3872)(-^ XcJl)K, we fix all the parameters 
obtained from the signal MC study after correcting the 
PDF shapes by applying MC/data calibration factors. 

To study background with a real J/il': we use large 
MC simulated B J/i/jX samples corresponding to 100 
times the integrated luminosity of the data. The non- 
J/ip (non-Xcj) background is studied using M« {Mj/.^^) 
sidebands in data. In the B — )• {xcJ7)K cases, the back- 
ground with a broad peaking structure is mostly due to 
the B ^ ^'(^ Xc]l)K, B ^ XcjK*, B ^ J/^jK* and 
B — ^ ip'K* decay modes. B — > Xc.'il)K produces 

peaks in both distributions (M-^^^-y and Mbc), while the 
other backgrounds are fiat in M^^^^ but peaked in Mbc. 
We determine the PDFs from the large MC sample. The 
fractions of the PDF components are floated in the flt, 
except for B — ^ ip'{-^ Xcil)K, whose fraction is con- 
trolled by fixing its ratio to the B — > V-''(^ Xc,}l)K sig- 
nal yield. For the combinatorial background, a threshold 
function, (Mp^^j^)^ x exp(a (Mp^^j.^ — Mth) + b (M^^^^ — 
Mth)2 + c (Mx^,^ -Mth)'), where Mth = 3.543 GeV/c^ 



(3.585 GeV/c^) for {M^^^^), is used for M^^^^^ and 

an ARGUS function [32[ is used for Mbc. The value of 
Mth is estimated from a MC study; its variation, which 
affects the signal yield in the fits, is incorporated in the 
systematic errors. The AE data sidebands are used to 
verify the background PDFs. The fractions for the signal 
and the background components are fioated in the fit. 

The results of the fits are presented in Figs. 1-3 and 
in Table HI The significance is estimated using the value 
of — 2 ln(£o/'Cmax) where >Cmax i^o) denotes the like- 
lihood value when the yield is allowed to vary (is set 
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FIG. 

cays; (a) M^^-^j distribution for Mbc > 5.27 GeV/c^, (b) 
Mbc distribution for 3.660 GeV/c^ < M-^^,^ < 3.708 GeV/c^ 
(-0' region, shown by red arrows), (c) Mbc distribution for 
3.805 GeV/c^ < M^^^j < 3.845 GeV/c^ (X(3823) region, 
shown by magenta arrows) and (d) Mbc distribution for 3.84 
GeV/c^ < M-^^,^ < 3.89 GeV/c^ (^■^'^2) region). The 
curves used in the fits are described in 1331. 



to zero). In the likelihood calculation, the x^ statistic 
uses the appropriate number of degrees of freedom (2 in 
case of B^ — > X{382i)K^ and 1 for the other decay 
modes). The systematic uncertainty, which is described 
below, is included in the significance calculation [s^l. We 
find a significant ip' signal in all considered channels. 
We also obtain evidence for the X(3823) in the channel 
B^ XcijK^ with a statistical significance of 3.8 stan- 
dard deviations (cr). The X(3872) signals are insignifi- 
cant. We estimate the branching fractions according to 
the formula B = 



^„ t' — ; here, Y is the yield, e is the 
reconstruction efficiency, is the secondary branching 
fractions taken from Ref. [9[ and Ngg is the number of 
the BB mesons in the data sample. Equal production of 
neutral and charged B meson pairs in the T(4S') decay is 
assumed. Measured branching fractions for the ip' £^re in 
agreement with the world average values for all the chan- 
nels d. We set 90% confidence level (C.L.) upper lim- 
its (U.L.) on the insignificant channels using frequentist 
methods based on an ensemble of pseudo-experiments. 

A correction for small differences in the signal detec- 
tion efficiency between MC and data has been applied for 
the lepton and kaon identification requirements. Uncer- 
tainties in these corrections are included in the systematic 
error. The 77 £+£' and D*+ -> D°{K-Tr+)Tr+ sam- 
ples are used to estimate the lepton identification correc- 
tion and the kaon (pion) identification correction, respec- 
tively. To estimate the correction and residual system- 
atic uncertainty for Kg reconstruction. — > D°(^' 
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FIG. 2: 2D UML fit projection for (Xc27)A'* decays: 

(a) M^^aT distribution for Afbc > 5.27 GeV/c^ (b) Mbc dis- 
tribution for 3.660 GeV/c^ < M^.^j < 3.708 GeV/c^ (i/-' re- 
gion), (c) Mbc distribution for 3.805 GeV/c^ < M^^^-y < 3.845 
GeV/c^ (X(3823) region) and (d) A'hc distribution for 3.84 
GeV/c^ < Mx^2^ < 3.89 GeV/c^ .^^i^'^^^ region). The 
curves used in the fits are described in [331 ■ 



KgTT^TT^)TT^ samples are used. The errors on the PDF 
shapes are obtained by varying all fixed parameters by 
±1(7 and taking the change in the yield as the system- 
atic uncertainty. The uncertainties due to the secondary 
branching fractions are also taken into account. The 
uncertainties of the tracking efiiciency and Ngg are es- 
timated to be 0.35% per track and 1.4%, respectively. 
The uncertainty on the photon identification is estimated 
to be 2.0%/photon. The systematic uncertainty asso- 
ciated with the difference of the 7r° veto between data 
and MC is estimated to be 1.2% from a study of the 
Xci(-> J/'il^l)K^ sample. 
We also perform a simultaneous fit to — !• 
{Xcil)K^ and B^ {Xcil)Kl, assuming that B{B^ 
X{382S)K^)/B{B° ^ X(3823)i^0) = B{B^ 
tjj'K^)/B{B° tP'K'^). The ijj' peak position and M^^,^ 
resolution arc common for both charged and neutral B 
candidates. From this fit, we estimate the significance for 
X(3823) to be 4.0 a (including systematic uncertainties). 
We determine the mass of the signal peak relative to the 
wcU-mcasurcd tJj' mass : 



M 



X(3823) 



-'"jf(3823) 



= 3823.1 ± 1.8 ± 0.7 MeV. 

Here, the first uncertainty is statistical and the second 
is systematic. Because of the mass-constrained fit to the 
Xci7 candidate, the systematic uncertainty of Mx(3823) is 
dominated by the additional photon's energy scale. This 
photon energy scale uncertainty is estimated by the dif- 
ference between the Xci J li^l candidates' mass with- 
out any constraint and the Xci nominal mass which 
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FIG. 3: 2D UML fit projection for B° ^ {Xc3l)K°s decays: 
(a) M^^i-y distribution for Mbc > 5.27 GeV/c^ (b) M^^^^ 
distribution for Afbc > 5.27 GeV/c^, (c) Mbc distribution 
for 3.660 GeV/c^ < M^^^^ < 3.708 GeV/c^ (V^' region), 
(d) Mbc distribution for 3.660 GeV/c^ < M^^^., < 3.708 
GeV/c^ {i>' region), (e) Mbc distribution for 3.805 GeV/c^ < 
Afx<:i7 < 3.845 GeV/c^ (X(3823) region), (f) Mbc distribu- 
tion for 3.805 GeV/c^ < M^^^^ < 3.845 GeV/c^ (X(3823) re- 
gion), (g) Mbc distribution for 3.84 GeV/c^ < M^^^j < 3.89 



GeV/c^ (X(3872) region) and (h) Mbc distribution for 3.84 
GeV/c^ < M-^^^^ < 3.89 GeV/c^ (X(3872) region). The 



curves used in the fits are described in [3S 



results in 0.7 MeV as the Afx(3823) systematic error. In 
order to estimate the width, we float this parameter and 
find no sensitivity with the available statistics: the width 
is 1.7 ± 5.5 McV. Using pseudo-experiments generated 
with different width hypotheses for the X(3823), the U.L. 
at 90% C.L. on its width is estimated to be 24 MeV. 

The mass of the X(3823) is near the potential model 
expectations for the l^Dj states' centroid mass: the Cor- 
nell [l^l and the Buchmiiller-Tye [2^ potentials give 3810 
MeV/c^. Other models predict the mass of -02 (the 
l^D2 cc state^^ having J^'^ = 2 — ) to be 3815-3840 
MeV/c^ di-ial. The X(3823) mass agrees quite well 
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FIG. 4: 2D UML fit projection of M-^^-^^-, distribution for tfie 
simultaneous fit of — ^ {Xci'y)K^ and — > (xci7)-K^s 
decays for Mbc > 5.27 GeV/c^. Tlie curves used in tlie fits 
are described in [3^ . 



TABLE I: Summary of tlie results. Signal yield (Y) from 
the fit, significance (5) with systematics included, corrected 
efficiency (e) and measured B. For B, the first (second) error is 
statistical (systematic). Here, in the neutral B decay case, the 
K% — >■ ■k'^-k~ branching fraction is included in the efficiency 
(e) but the factor of 2 (for K'^ — >■ K% or K^) is separately 
taken into account. 



Decay 


Yield (F) 5(o-) 


.(%) 


Branching fraction 




XcJ7)^* 




s(io-'') 




193.2 ±19.2 14.8 


8.6 


7.7 ±0.8 ±0.9 




59.1 ±8.4 7.8 


6.0 


6.3 ±0.9 ±0.6 












50.3 ±7.3 7.2 


5.1 


6.8 ±1.0 ±0.7 




12.9 ±4.4 2.9 


3.5 


4.7 ±1.6 ±0.8 




X(3823)(^ Xcn)K^ 


e(io-'') 




33.2 ±9.7 3.8 


10.9 


9.7 ±2.8 ±1.1 




0.3 ±3.9 0.1 


8.8 


< 3.6 




X(3823)(^ XcJ7)^" 






Xcl 


3.9 ±3.4 1.2 


6.0 


< 9.9 


Xc2 


5.3 ±2.9 2.4 


5.0 


< 22.8 




X(3872)(-> Xc.i7)A'* 




Xcl 


-0.9 ±5.1 


11.1 


< 1.9 


Xc2 


4.7 ±4.4 1.3 


9.3 


< 6.7 




X(3872)(^ XcJ7)^" 






Xcl 


4.6 ±3.0 1.6 


6.2 


< 9.6 


Xc2 


2.3 ±2.2 1.1 


5.2 


< 12.2 



with these models. In addition, since no peak has been 
seen around X(3823) in the DD final state [H], one ex- 
pects that ip2 does not decay to DD [s^. The ratio 
Rb = < 0-41 (at 90% C.L.) is consistent 

with the expectation {Rb ^ 0.2) for -02 [391 |40|. The 
limited data preclude an analysis of the decay product 
angular distribution to determine the J^*^ assignment of 
the X(3823). The product of branching fractions for the 



X(3823) is approximately two orders of magnitude lower 
than for the ip' , as shown in Table HI it is consistent with 
the interpretation of the X(3823) as "02(1^-02), whose 
production rate is suppressed by the factorization j4l| in 
the two-body B meson decays. 

In summary, we obtain the first evidence of a nar- 
row state, X(3823), that decays to Xci7 with a mass 
of 3823.1 ± 1.8(stat) ± 0.7(syst) MeV/c^ and a sig- 
nificance of 3.8 cr, including systematic uncertainties. 
We measure the branching fraction product B{B^ 
X{3823)K^)B{X{3823) xal) = (9.7 ± 2.8 ± 1.1) x 
10~^. No evidence is found for X(3823) -> Xc27 and 
we set an U.L. on its branching fraction product B as 
well as the ratio Rb = < 0.41 at 90% 

C.L. The properties of the X(3823) are consistent with 
those expected for the tp2 {^^D2 cc) state. We also deter- 
mine the U.L. on the product of the branching fractions, 
B{B^ X{3872)K^)B{X{3872) xal) < 1-9 x IQ-^ 
at 90% C.L.; this is less than one quarter of the corre- 
sponding value in X(3872) J/f/jTr+Tr" 0. In addition, 
we obtain constraints for all of the Xc]jK final states in 
neutral and charged B decays. Our results show that the 
X(3872)'s C-odd partner production in the two-body B 
decays or its decay into XcJ7 are found to be considerably 
suppressed. 
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